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Abstract — In this paper, we study network coding capacity 
for random wireless networks. Previous work on network 
coding capacity for wired and wireless networks have focused 
on the case where the capacities of links in the network are 
independent. In this paper, we consider a more realistic model, 
where wireless networks are modeled by random geometric 
graphs with interference and noise. In this model, the capacities 
of links are not independent. We consider two scenarios, single 
source multiple destinations and multiple sources multiple 
destinations. In the first scenario, employing coupling and 
martingale methods, we show that the network coding capacity 
for random wireless networks still exhibits a concentration 
behavior around the mean value of the minimum cut under 
some mild conditions. Furthermore, we establish upper and 
lower bounds on the network coding capacity for dependent 
and independent nodes. In the second one, we also show that the 
network coding capacity still follows a concentration behavior. 
Our simulation results confirm our theoretical predictions. 

I. Introduction 

Network coding was originally proposed by Ahlswede et 
al. in [1]. Unlike traditional store-and-forward routing al- 
gorithms, in network coding schemes, intermediate nodes 
encode their received messages and forward the coded mes- 
sages to their next-hop neighbors. It has been shown that 
network coding can improve the network capacity, even by 
using simple linear or random codes [8], [9], [11], [12]. 
In most studies of network coding, network topologies are 
assumed to be known. 

In [17], [18], the authors studied network coding capacity 
for weighted random graphs and random geometric graphs. 
In the random graph model, each pair of nodes are connected 
by a bidirectional link with probability p < 1 indepen- 
dently [4], [10]. The capacity of each link is assumed to 
be i.i.d. according to some probability distribution. In the 
random geometric graph model, two nodes are connected to 
each other by a bidirectional link only when their distance 
is less than a predefined positive value r, the characteristic 
radius [15]. Each link has a unit capacity. For these two 
types of random networks, the authors showed that the 
network coding capacity is concentrated at the (weighted) 
mean degree of the graph, i.e., the (weighted) mean number 
of neighbors of each node. Essentially, the results reveal 
a concentration behavior of the size of the minimum cut 
between two nodes in random graphs or random geometric 



graphs. Similar problems have been studied in the literature, 
e.g., [6] and references there. In [3], the authors studied a 
generalized random geometric graph model, where two nodes 
are connected by a bidirectional link with probability 1 if 
their distance d is less than ro > and with probability 
p < 1 if r*o < d < r\. They obtained similar concentration 
results there. 

The geometric models in [3], [17], [18] assume that a 
link exists (possibly with a probability) between two nodes 
when the nodes are within each other's transmission range. 
Although each link has a direction, as all links are bidirec- 
tional (i.e., the link implies the existence of the link 
(j, i)), the model in fact leads to an undirected graph and 
considerably simplifies the resulting analysis. In addition, 
interferences among wireless terminals were not considered 
in [3], [17], [18]. Nevertheless, in wireless networks, due to 
noise, interference, and heterogeneity of transmission power, 
significantly more sophisticated models for link connectivity 
are needed. For instance, a widely-used model for wireless 
communication channels is the Signal-to-Interference-plus- 
Noise-Ratio (SINR) model [16], [19]. In this paper, we study 
the capacity, i.e., the size of the minimum cut, of random 
wireless networks under the SINR model. 

Since how to apply the network coding with noisy links 
is still an open problem, we assume that as long as the 
SINR of a link (3ij is greater than or equal to a 

predefined threshold (3, then node i can transmit data at rate 
R packets/sec to node j without any error. That is links are 
noise-free once the SINR condition is met. In other words, 
we view the network coding as operation on a higher layer 
in the network communication stack, and assume there is 
an error correcting code at the lower layer which corrects 
errors on the links once the SINR threshold is met. Then, in 
this model, each link is indeed directional (not necessarily 
bidirectional), and the capacities of different links are not 
independent. We will show that the capacity still has a sharp 
concentration when the scale of the network is large enough. 

This paper is organized as follows. In Section II, we 
describe the random wireless network model. In Section 
III, we study the network coding capacity for a single 
source and multiple destinations transmissions. Specifically, 
we investigate two cases. In the first one, all nodes have 



the same transmission power, and in the second one, the 
transmission powers are heterogeneous. We use different 
techniques for these two cases and show that the network 
coding capacity has a concentration behavior in both cases. 
In Section IV, we extend our result to multiple sources and 
multiple destinations transmission problem. In Section V, we 
present some simulation results, and finally, we conclude this 
paper in Section VI. 

II. Random Wireless Networks Model 

We use the following model for random wireless networks. 
Assume 

(i) X = {Xi, X2, X„} is a set of i.i.d. two-dimensional 
random variables according to a homogeneous Poisson 
point process in the two-dimensional unit torus, where 
Xi denotes the random location of node i, and n is the 
total number of nodes. 

(ii) Each node i has a transmission power Pi, which follows 
a probability distribution fp(p), p G \p m m,Pmax], 
where < p min < p max < 00. 

Here, the existence of a link from node i to node j depends 
on the ability to decode the transmitted signal from i to j, 
which is determined by the Signal-to-interference-plus-noise- 
ratio (SINR) given by 



P l L{d lJ ) 



(1) 



where Pi is the transmission power of node i, dij is the 
distance between nodes i and j, and TVo is the power of 
background noise. The parameter 7 is the inverse of system 
processing gain. It is equal to 1 in a narrow-band system and 
smaller than 1 in a broadband (e.g., CDMA) system. The 
signal attenuation function L(-) is a function of the distance 
d,^ = ||Xj — Xj||, where || • j| is the Euclidean norm, and 
is usually given by L(dij) = cd~j° for some constants c and 
2 < a < 4. 

Under the SINR model, the transmitted signal of node 
i can be decoded at j if and only if /3y > /?, where (3 
is some threshold for decoding. In this case, a link 
is said to exist from i to j. Note that even if (3ij > (3, 
(3ji > (3 may not hold and thus the link (j, i) may not exist. 
Thus, the graph resulting from the SINR model is in general 
directed. It is clear that link is bidirectional if and only 
if min{Pij , Pji} > (3. Denote by G(X ,V ,7) the ensemble 
of random wireless networks induced by the above physical 
model, where V = {Pi, P2, P n } represents the set of 
transmission power. 

For transmission power P and signal attenuation function 
L(-), we assume 

(i) Pmin > (3N ; 

(ii) Pr(P = Pmm ) > 0, Pr(P = Pmax ) > 0, 

(iii) L(x) is continuous and strictly decreasing in x 

for technical and practical reasons. In the remainder of this 
paper, under different circumstances, we may add further 
constrains on Pi. 




Fig. 1. 
graphs 



Single-source single-destination transmission in directed SINR 



The sum £ fc?y L(d k j) = J2 k ^ L (W X k - X,-||) is a 
random variable depending on the locations of all nodes in 
the network. Define 



for all j, 



for all j. 



(2) 



(3) 



To study the asymptotic network capacity, we will let the 
number of nodes n go to infinity. Since the region is fixed, 
this corresponds to a dense network model [?], [15]. Another 
widely used model is the extended network model [7], [13], 
in which the number of nodes and the area of the region both 
go to infinity while the ratio between them — the density of 
the network, is kept as a constant. Both models are widely 
used in the literature. We will focus on the former one in 
this paper. 

III. Network Coding Capacity for Single Source 
Transmission 

A. Capacity of a Cut 

Let Cij be the capacity of a link We will specify 

the value of Cij later for different scenarios. Consider 
a single-source multiple-destination transmission problem. 
Let s be the source node. Suppose there are I destination 
nodes, t\,...ti, and m relay nodes, u\,...u m . Denote the 
set of the destination nodes and relay nodes by T and 1Z, 
respectively. Fig. 1 illustrates an example of single-source 
single-destination transmission. 

Let the capacity of the link from the source s to each relay 
node Ui be C S i, i = 1, ...,m, the capacity from relay node 
u.i to another relay node uj be Cy, i 7^ j,i = 1, m, j = 
l,...,m, and the capacity from each relay node it, to each 
destination node tj be Cu-, i = l,...,m, j = Unlike 
random geometric graph models studied in [3], [17], [18], the 
capacities in our model are not symmetric nor independent 
in general. 

Since in our random SINR wireless network model, there 
are two sources of randomness: one is the random location 
of each node and the other is the random transmission power 




Fig. 2. An s-i-cut for the single-source single-destination transmission in 
directed SINR graphs 



of each node. We use Ex and Ep to denote the expectation 
operation with respect to each probability measure respec- 
tively. 

Let C be the expected capacity of a link which is 
defined as 



C 



ExEp[Ci 



CijdF 0ij (r), 



(4) 



where F^A-) is the c.d.f. of /3y, which is determined by 
/p(-), the distribution of X, and path-loss function L(-). 

Now define an s-t-cut of size k for a pair of given source 
s and destination t G T as a partition of the relay nodes 



into two sets 14 and V k c , such that |Vi-| = k, \V k c 



k, 



V k U14 C = K and V k fl V k c 
is shown in Fig. 2. Let 



An example of an s-t-cut 



c k = J2 c « + E E c * + E cj*. © 

UjEV k 



UiEV k UiEV£ 



then C k is the capacity of the corresponding s-t-cut. Al- 
though Cfc is a sum of dependent random variables, we still 
have 



in the SINR model, Ck is a sum of dependent random link 
capacities. Instead, we employ coupling, martingale methods 
and Azuma's inequality [14], [15] to solve the problem for 
different cases. 

Note that when 7 = 0, i.e., there is no interference in the 
networks, the capacities C S i for i = k+1, m are mutually 
independent, as well as the capacities dj for any fixed 
i = 1, k with j = k + 1, m or t. In this case, although 
the link capacities are still asymmetric, J2u gv anc ^ 
J2u EVu{t} Cji for j <G Vk become sums of independent 
random variables. Thus we can apply methods similar to 
those used in [3], [17], [18] to obtain the same concentration 
results. 

B. Constant Transmission Power 

Consider the scenario when all nodes transmit with a 
constant power Pq and denote the model by G(X , Po,7)- 
The SINR of link (i, j) in this case, flij, can be rewritten as 

L(dij) 



fl, 



No/Po+lEk^Hd k] ) 
L(dij) 



(7) 



N /P Q + 1 J(j)- 1 L{d ij) 

Assume when 0^ > (3, the link (i, j) has a capacity R, 
i.e., node i can transmit data at rate R packets/sec to node 
j without any error. Then, we can define as 

R Pi] > P, ( o\ 
o Pa < p. w 

Note that when the wireless channel is Gaussian channel, 
the capacity of link (i, j) is [5] 

1 



Ci 



c, 



log(l + /%) Pij>P, 

Pa < P- 



(9) 



E[Ck] 



Our results in this subsection do not rely on any particular 
expression of Cij, and thus they hold for Cy defined by (O 
as well as for dj defined by (0. Nevertheless, since we 
consider the application of network coding, it would be more 
appropriate to focus on the R — capacity ([S]), rather than 
the capacity of Gaussian channel. 



— Ex Ep [Ck] 

= ^ E x Ep[C si ] + E E E xE P [C jt ] jyy Because^of the i.i.d. distribution of X./s, given X^, 



Note that P^ and thus C tJ 



are determined by L(dij) and 



tnev; 
+ J2 E xE P [C lt ] 

[to + k(m - k)]C, 



djj 's are independent for all i ^ j. Given node j, let 



E[L] ± E Xi [L(dij)], 



(10) 



(6) then 



and consequently E[Ck] = E[C m -k] for k = 0,1,..., m, 
and E[C ] < E[d] <■■■< E[C {m/2] }. 

To show the capacity of any source-destination pair con- 
centrates at some value, we will first show that for such a 
source-destination pair, the capacity of any s-t-cut of size k 
concentrates at its mean value. Similar results were proved 
in [3], [17], [18], where the capacities of the links that 
originate from the same node are i.i.d. Nevertheless, the 
methods used in [3], [17], [18] do not apply here, since 



E[J(j)} = E 



Y,L{d, 



= (n-l)E[L]±E[J]. (11) 



Since our model is a dense network model and the area 
of the region is fixed, E[L] = E[L(dij)] is a constant and 
E[J] = (n — 1)E[L] scales with n. For different j's, it is 
clear that J(j)'s are not independent, however, they have 
the same sharp concentration behavior in large scale wireless 
networks. This is established in the following lemma. 



Lemma 1: Suppose there are n nodes in the network, then 

Pr(J(j) < (1 - ei)E[J\) = O (±) , (12) 

Pv(J(j)>(l + e[)E[J]) = o(-^j, (13) 
for all j = 1,2, ...,7i, where e\ = J 7~^§pi an d e i = 



and 



6 In n 



(n-l)-E[ip 

Proof: Given any node j, because J(j) = 
and L(dij) are i.i.d. for all i ^ j, by the Chernoff bound [2], 
[14], we have 

E[J]el 



Pr(J(j) < (l-ex)E[J]) < exp 
= exp 



(n-l)E[L]el) 



.14) 



and 



Pr(J(j)>(l + ci)^W) <exp 
= exp 



(n - 



15) 



Substituting £1 = yjjz*&jht and ej = ^ into 
(fl4l and ([TBI l. we obtain ( fT2l and d 13t . respectively. □ 

Lemma [T] shows that when the network is large, i.e., n is 
sufficiently large, the interference at each node concentrates 
at 7(n - l)E[L] = Q(1)E[L]. The reason for this is the 
uniformly (asymptotically Poisson) random distribution of 
the nodes. 

Now define two other types of SINR models G'(X, Pq, 7) 
and G" (X 1 Pq, 7) which are coupled with G(X, Po, 7) such 
that they have the same point process X and constant 
power P - Let the SINR of link in G'(X,Po,j) and 
G"(Af,P ,7) be 

L(dij) 



and 



N /P + (l + e[) 1 E[J]- 7 L(d ij ) 

Ljdij) 

N Q /P Q + {l-e 1 ) 1 E[J]- 1 L{d ij y 



(16) 
(17) 



respectively. 

Let and C^ be the capacity of link in 

G'(X,P ,j) and G"{X,P ,j), respectively. Since e x -> 
and e'j — * as n — > 00, Cy and C 4 " are asymptotically equal 
to Cij. 

The following lemma establishes a concentration result for 
Ck with constant transmission power by coupling methods. 

Lemma 2: For any < e < 1, the capacity of an s-t-cvX 
of size k, k = 0, 1, ...,m, satisfies 



Pr(C fc < (l-e)E[C' k }) < exp 



l-O 



where S[C1] = [m + k(m — k)]C and C is the average 
link capacity in G'(X,Pq,j), and 



Pr(C fc > (l+e)E[(%\) < exp - 



E[C'{\ 




l-O 



where -E^C^'] = [m + fc(m — fc)]C" and C" is the average 
link capacity in G"(X, Po,7). 

Proo/> Since for all j, {J(j) > (1 - ei)^[J]} and 
{^(j) < (1 + e 'i)E[J]} are both increasing events^ By the 
FKG inequality [2], [13], [15], we have 



Pr f]{J(j) > (1 - >I1 Pr ( J W ^ ( X - 

Vj=i / 3=1 

l-Ol 



= 1-01 - 

n 



where the first equality is due to Lemma [T] and 

/ n \ n 

Pr f| { J(j) < (1 + el)S[J]} >II Pr(J(7) < (1 + £ 'i)£[J] 

\j = l / 3=1 

l-Ol 



=1-01 - 



This implies that Cy is stochastically lower bounded by 
and stochastically upper bounded by C," with probability 
1 — 0(— ). Hence, in order to show (fT8l ) and (fT9l) , it suffices 
to show 



Pr(C' k <(l-e)E[C k })<exp 



and 



Pr(C£ > (1 + e)P[C£]) < exp - 



En 



(18) 



(20) 



(21) 



In G'(X, P , 7) and G"(X, P , 7), the SINR of link (i, j) 
is given by ([Tol l and (fTTT i, respectively, and because d^'s 
for a given i are independent, by applying the Chernoff 
bounds, we obtain (l20b and (|2TT i. □ 

Since and C," are asymptotically equal to Cij, E[C' k ] 
and -E[C^,'] are asymptotically equal to E[C k ]. Consequently, 
Lemma 12 shows that C k concentrates at E\C k ] asymptoti- 
cally almost surely (a.a.s.). 

Now, let C S: t be the minimum cut capacity among all s- 
t-cuts, i.e., 

C s .t = min Cfc. (22) 

0<fe<m 

'in context of graph theory, an event A is called increasing if Ia(G) < 
Ia(G') whenever graph G is a subgraph of G' , where I a is the indicator 
function of A. An event A is called decreasing if A c is increasing. For 
details, please see [2], [13], [15]. 



For the given source node s and the sets of destination nodes 
T = {t±,...,ti} and relay nodes 1Z = {ui, u m }, define 
the network coding capacity as 



C s ,t = min C s .t- 
ter 



(23) 



That is because for one source and multiple destinations, the 
capacity of network coding depends on the minimum cut 
among all the destinations. 

In the following, we show that when the number of relay 
nodes m is sufficiently large, the network coding capacity 
C S .T concentrates at E[Cq] = mC with high probability. 

Theorem 3: When n is sufficiently large, with high prob- 
ability, the network coding capacity C s .r satisfies 

Pr(a. r > (1 - e' a )E[C Q ]) = 1-0 (-L\ , (24) 



where e' 



2a lr 



for a > and E[C Q } = mC. 



Proof: Since the C^'s are asymptotically equal to C^-'s, 
in order to show d24l i. it is equivalent to show 

Pr(C 8 , r >(l-«]) = l-0 

Since E[C' k ] > E[C ] for any k = 1, m, 
Pr(C s , 4 < (1 - e' a )E[C' Q ]) < Pr(C M < (1 - e' a )E[C k ,]), 

for any t <G T, where fc' is the size of the minimum s-i-cut. 
By dT8b of Lemma |2j we have 



Pr(C s , t < (l-e' a )E[C' k ,])<exp 



<exp 



e' Q 2 [m + fc'(m-fc')]C" 



2a lr 



E[C _ 



-, since C" and C are asymptot- 



By choosing e' a - 

ically equal, we have for any t E T. 

Pr(C7 Sil <(l-eMg) = o(l 
By the union bound, we have 



Pr(C s , T < (1 - 4)^[^])<X; P 1 ^.* ^ (! " 



ter 



=0 



□ 



Theorem 4: When n is sufficiently large, with high prob- 
ability, the network coding capacity C s .r satisfies 

Pr(C s , T < (1 + %)E[C ]) = 1-0 (±\ , (25) 



where e 1 ' 



3a In tn 
E[C ] 



for a > and S[C ] = mC. 



Proof: Since the C^'s are asymptotically equal to C^'s, 
in order to show d25l >. it is equivalent to show 

Pr(C s , T <(l + eMa = l-o(i 

To show this, it is sufficient to consider a particular cut for 
a pair of the source and one destination, e.g., an s-i-cut 
separating the source s from all the other nodes. 

Pr(C s ,r > (1 + C)E[C'o]) < Pr {C s ,t > (1 + &)E[(%]) 

< Pr [Y,C si >{l + e^)E[CZ] 

\i = l 

= Pr(C > {1 + <&)E[CS]) 
e'fmC" ' 

< exp 



= O 



where the last inequality follows from (fT~9b of Lemma [2] 

□ 



C. Heterogeneous Transmission Powers 

In this subsection, we consider the case where the trans- 
mission power of each node is random rather than a constant, 
but the capacity of a link is a constant R, which is 

independent of the SINR fyj, when /3ij > (3. In this case, 
Pij can be rewritten as 

PiL{dij) 



ft 



P t L(d tJ ) 



(26) 



Because Pi's and X,'s are both i.i.d., using the same 
method, we can prove the following lemma: 

Lemma 5: Suppose there are n nodes in the network, then 



Pr(7(i) < (1 - e 2 )E[I}) = O I i 



and 



Pv(I(j)>(l + e' 2 )E[I}) = -j 



(27) 



(28) 



for all j = 1,2, 7i, where 62 



4 In 71 „ nr i / 

(n-l)S[P]B[L] anQ £ 2 



6 In n 



(n-l)E[P]E[L] ■ 

Even though we have concentration results for I (J), we 
cannot employ the same coupling methods as in the previous 
section. This is because in G'(X, Po,j) (or G"(X, Pq, 7)), 
the C^-'s (respectively, C^-'s) are independent for all j 7^ i 
for given i. In our new case, however, this independence does 
not hold because all Cj/s depend on transmission power Pj. 
To deal with this dependence, we use martingale methods 
and Azuma's inequality to solve our problem. 

Theorem 6 (Azuma's Inequality [2]): Let Z$, Z\, be a 
martingale sequence such that for each i = 1,2,...,, 

\Z t - Zi-\\ < Ci 



almost surely, where c, may depend on i. Then for all n > 
and any A > 0, 



Pr(Z n > Z Q + A) < cxp <{ - 77T ^ ^ j> , (29) 



and 



Pr(Z n < Z - A) < cxp < - 




Proof: Please see e.g. [2]. 



(30) 

□ 



To use Azuma's inequality, we need to construct a mar- 
tingale. A common approach to obtain a martingale from a 
sequence of random variables (not necessarily independent) 
is to construct a Doob sequence. More precisely, suppose we 
have a sequence of random variables Yi, Y 2 , Y n , which 
are not necessarily independent. Let S = J^™ Yi and define 
a new sequence of random variables {Zi : i = 0, 1, ...,n} 
by: 

Za = E[S] 

Zi = E Yz+u ... tYn [S\Y u ...,Y^, i=l,2,...,n. 

(31) 

Then {Zi : i = 0, 1, n} is a martingale and Z n = S. 

If we are able to upper bound the difference \Zi — Zi-\ | for 
all i by some constant, then we can apply Azuma's inequality 
to obtain some bound on a tail probability. For example, if 
Yi's are independent, a simple upper bound for \Z± — Zj-il 
is any upper bound on |Y^|. However, as long as the Yi's are 
dependent, which is the case in our model, we cannot bound 
\Zi — Zi-\\ in this way. In this case, we need to understand 
the properties of the Yi's to see if we can bound \Zi — Zi-\\. 
We approach our problem by following this idea and using 
the next corollary. 

Lemma 7: For n > 1, given a sequence of random vari- 
ables Y\, Y 2 , Y n , which are not necessarily independent, 
let S = Yn^i- If f° r an y ViiVi S Di, where Di is the 
support of Yi, 

\E[S\Y 1 ,...,Yi_ 1 ,Y i =y i }-E[S\Y 1 ,...,Y i _ 1 ,Y i = y(}\<c i , 
almost surely, where may depend on i, then for any A > 0, 

and 



Pr(S>E[S] + X)<eiq ? \-^ w 



(32) 



(33) 



Proof: We prove this corollary for the case of discrete 
random variables. For continuous random variables, the proof 
is similar. 

Define a Doob sequence with respect to {Yi : i = 
1, ...,n} as in QTl ). To simplify the notation, we will write 

E Yi+1 ,.,.,Y n [S\Y 1 ,...,Y i \ as E[S\Y U ...,Yi\ when there is 
ambiguity. 

By the total conditional probability theorem, we have 

Zi^E^Yu.^Yi-i] 



= Y E[S\Y 1 ,...,Y i - 1 ,Yi = y]Pr(Yi = y\Y 1 ,...,Y i - 1 
yea, 

and 

Z, = EiSlYt, ...,Yi\ 

= YI E[S\Y u ...,Y i \Pv(Yi = y\Y 1 ,...,Y i _ 1 ). 

yeDi 

Therefore, 

\Zi-Zi_ x \ 
= \E[S\Y 1 ,...,Y}-E[S\Y 1 ,...,Y_ 1 }\ 

= i Y E l s \ Y ^ Pr (^ = y\ Y ^ 

- Y E[S\Y 1 ,...,Y i - 1 ,Y i = y}Pv(Y i = y\Y 1 ,...,Y i . 

< Y \E[S\Y 1 ,...,Y i ]-E[S\Y 1 ,...,Y^ lt Y i = y]\ 

yeDi 

■Pr(Y i =y\Y u ...,Yi_ 1 ) 

yeDi 

= Ci. 

Since {Zi : i = 0, 1, n} is a martingale with bounded 
difference of \Zi — Zi_\\, we can apply Azuma's inequality 
to obtain d32]i and £2}. □ 

Now consider G'(X, V, 7) and G"(X, V, 7) coupled with 
G(X, P, 7) such that they have the same point process X 
and powers V. Then, the SINR of link in G'(X,V,j) 
and G"(X,V,j) are 



and 



Pij 
0"- = 



11 N a + {l + e' 2 ) 1 E[I]- 1 P l L(d l3 ) 
PiHdij) 



v N Q + {\~e 2 ) 1 E[I]~ 1 P l L{d l3 y 



(34) 
(35) 



respectively. 

Let C'^ and C," be the capacity of link in 
G'(X,V,~/) and G'\X,V,i), respectively. Then, and 
Cj" are asymptotically equal to Cy. 

Assume that there exist r' mm > 0, r' max > 0, r^ m > 



and r" lax > as the solutions for 



and 



PminL(r m i n ) 


No + 7(1 - 


¥e' 2 )E[I] - lPmi n L{r' min ) 




PmaxL(r max ) 


iV + 7 (l- 


^ e ' 2 )E[I] - lPmaxL{r' max ) 




PminL{r min ) 


N + 7(1 - 


- e 2 )E[I] - JPminLiCJ 




PmaxL(r max ) 



N + 7(1 - e 2 )E[I] - lp max L{r' 7 ' nax ) 
respectively. That is 

-if P N 0+1 (l + e' 2 )E[I] 



r' = U 



1 + 7/3 



= L~ 





/3 


N + j(l + e( 2 )E[I] 


1 


+ 7/3 


Prnax 




P 


N Q + j(l - e 2 )E[I] 


1 


+ 7/3 


Pmin 




P 


No + i(l - e 2 )E[I] 


1 


+ 7/3 


Pmax 



r' ■ 

mm' 



Since L(-) is continuous and strictly decreasing, r 
r'max> Cm and r 'max are a11 unique. In G'(X,V,i) 
(G"(X, V, 7)), any node inside the circle centered at with 
radius r' min (r" lin ) is connected to node i by a bidirectional 
link; while any node outside the circle centered at X^ with 
radius r' max (r^ ax ) is not connected to node i. 



Let A(X-i,r min ,r ma 



and inner radius 



Denote by N(r' r 
of nodes in „4(X. 



max-i 

^ and A{yi u r'^ m ,r'^ ax ) be the 
/ 

min 

mi n and outer radius r" lax , respectively. 

) the number 



max' 



mm ' maa; 



and N{r' 7 ! r 



mzn ' max 



1 1 ^min ' ' max 



and A(Xi 



r. 



,.11 



r' ) and N(r" r" ) 
min ' max J aIIU V min ' max ) 



spectively. It is clear that N(r' Tl 
have Poisson distribution with mean nirfr!^ 
™ {r'max ~ Cm), respectively. 

Now suppose the signal attenuation function L(x) 
for some constants c > and 2 < a < 4. Then, 



J and 



and 



r' 2 

max 



111 



r' 2 - 

in/n 



7/2 



c(l 


+ lP)Pmin 


/3[iV + 


7(1 + 4)^1/]] 


c(l 


+ lP)Pmax 


/3[iVo + 


7 (1 + e 2 )£[/]] 


c(l 


+ lP)Pmin 


/3[iVo + 


7(1 - e a )£;[i]] 


c(l 


+ lP)Pmax 


/3[iV + 


7(1 - e 2 )£[/]] 




{.Pmin 1 Pmax 



V 



[Wo + 7(1 



e' 2 )£[/]] 2 « 



mm 1 Fmax ) 



where 



= (p 



2a 
max 



[No+j(l-e 2 )E[I]} 2a 
c(l+ 7 /3)" r 



p 2a 

Jr'rmn 



P 



^From d36T l and (137) . we can see that both ^^(r'^ 



and mr(r 



J scale with n as 11,-1 , since 



maa: mm / 

scales linear with n. Now assume that there exists a constant 



i] > independent of n such that 

N {r'min, r ' m ax) <V> and ^0" 



mm i max 



) < V 09) 



hold a.a.s. This assumption actually puts a constraint on the 
transmission power since it needs to scale (if it scales) with 
n so that (|39l is satisfied. For example, we may choose 
r/ = 1 and the transmission power P scales with n so that 
B{pmimPmax) = ©(t^tt)- Note that r' min and r' max are 
asymptotically equal to Cin and r" n ax , respectively. 

The following lemma establishes a concentration result 
for Ck with heterogeneous transmission power and constant 
capacity by coupling methods and Azuma's inequality. 



Lemma 8: For any < e < 1, when n is sufficiently large 
and (|39l is guaranteed, with high probability, the capacity of 
an s-i-cut of size k, k = 0, 1, m, satisfies 



Pr(C fc < (l-e)£[C£])<expj- 



m + k(m - k))C' 
2{rj + l) 2 R 2 



(40) 

where E[C' k ] = [m + k(m — k)]C and C is the average 
link capacity in G '{X ,V ,7), and 



Pr(C fc > (l+e)£[C£]) < exp - 



[m + k(m - k)]C" 

2( V +l) 2 R 2 

(41) 

where -E^C^'] = [m + fc(m — fc)]C" and C" is the average 
link capacity in 7>, 7). 

Proof: By Lemma [5J for all j, (1 - e 2 )E[I] < < 
(1 + e' 2 )E[I] holds a.a.s. It is clear that CV, is stochastically 
lower bounded by C'^ and stochastically upper bounded by 
C"j almost surely. Hence, in order to show d40t and i41i . it 
suffices to show 

Pr(^<(l-e)i«)<exp' \rn + k { m - k))G^ 



and 



Pr(£7£ > (He)£[C[]) < exp 



2(77+ l) 2 i? 2 

(42) 

[m + k(m - k)]C" 2 e 2 
2(77+ l) 2 i? 2 



(43) 



To show (02), we use martingale methods. Let Yi 



s(fe+l) 
l(fe+l) 

and Y, 



s(k+2) ' ■■ 

^l(fe+l) ' ^m-fc+2 = Ci(fe+2) i •"! ^m-fc+fe(m-fc) 



'-'km' 



C 2 n Ym — k+k(r 



m-k+k(m~k) + l — C[ t , l / m -fc+fc(m-fc)+2 — 

= C^. t . Define a Doob sequence 
{Zj : i = 0, ...,m + fc(m — fc)} with respect to 

{Yi, i = 1, 2, m + fc(m — fc)} as 

(36) j Z = E[C' k ] 

\ Zi = ElC'^Yt, ...,Yi], i = 1,2,..., m + k{m- k). 

(37) Then {Z t : i = 0, ...,m + fe(m — k)} is a martingale and 

Zm+k{m-k) = Cfc- 

Since when j^u and j 7^ v, C-j is independent of C' uv , 
we have only the dependence among C^'s for all j ^ i 

(38) w j tn gi ven ^ However, the distance d^'s are independent 
for all j 7^ i with given i. When dij < r' min , C-j = R, and 

) when > r' „„, C" = 0. Moreover, the number of nodes 



within the annulus *4(X i5 r' min , r' max ) is upper bounded by 
the constant r/ a.a.s. Therefore, we have 

\E[C' k \Y x , ...,Yi_ lt Yi = Vi ] - E[C' k \Y x , ...,^-1,^ = y[]\ 

< (»7+l)-R 

a.a.s., where yi and are either or R. Applying the result 
of Lemma [7] we have ( 142) . In the same manner, we can 
show that (03) holds. □ 

In the following, we show that as the number of relay 
nodes m is sufficiently large, the network coding capacity 
C s .r concentrates at E[Co] — mC with high probability. 
The proofs are similar to those for Theorem [3] and Theorem 

El 



Theorem 9: When n is sufficiently large, with high prob 
ability, the network coding capacity C s .r satisfies 

I 



Pi(C. T >(l-e a )E[C ]) = l-O 



(44) 



where e a = ^rg~r V 2am In to for a > and -B[Co] = mC. 

Proof: Since C^'s are asymptotically equal to Cy's, in 
order to show d44l i. it is equivalent to show 

Px(C SiT >(l-e a )E[C^}) = l-o(-L 

Since E[C' k ] > E[C ] for any k = 1, to, 
Pr(C M < (1 - e a )£[C ]) < Pr(C M < (1 - e a )£[C£,]), 

for any t € T, where fc' is the size of the minimum s-t-cut. 
By d40b of Lemma [8j we have 

Pr(C s , t < (1 - e a )E[C k ,]) < exp <j 2(r? + l) 2 i? 2 

e 2 a mC' 2 
< exp < 



2(r) + l) 2 B? 



By choosing e a = ^ft^ f V 2aTO In to for a > 0, since C" 



and C are asymptotically equal, for any t £ T, 
Pr(C s , t <(l-e a ) J E[C^]) = of^ 
By the union bound, we have 

Pr(C s>r < (1 - e a )E[C' ]) < ^Pr(C M < (1 - e„)£[CS]) 



/ 



□ 



Theorem 10: When n is sufficiently large, with high prob- 
ability, the network coding capacity C s .r satisfies 



Pr(C s ,T < (l + e a )E[C Q }) = 1-0 



1 



(45) 



where e a = ^gjg^ V 2«to In to for a > and -B[Co] = mC. 

Proof: Since C^'s are asymptotically equal to Cy's, in 
order to show d45l l. it is equivalent to show 



Pr(C i ,r<(l + ea)«]) = l-0 — 

y TO 

To show this, it is sufficient to consider a particular cut for 
a pair of the source and one destination, for instance, an 
s-t-cut separating the source s from all the other nodes. 

Pt(C s ,t > (1 + e a )E[C'o]) < Pr (C s ,t > (1 + e a )E[C'o}) 

/ m 

< Pr [Y^Csi > (l + ec)E[CS] 
\i=i 

= Pr(C > (l + e a )E[CS]) 
eimC" 2 



< exp 



2(77 + l) 2 i? 2 



where the last inequality follows from (T4Tb of Lemma [8] □ 



IV. Network Coding Capacity for 
Multiple- Source Transmission 

In this section, we study network coding capacity for 
multiple sources and multiple destinations transmission. We 
assume the same notation as in Section III. However, instead 
of having a single source, we have s > 2 sources. Denote 
by S = {s\, . . . , Sh} the set of source nodes. Assume there 
is no correlation among the set of sources S. Now we can 
define an S-t-cut of size k between the set of sources S and 
one destination t E T as a partition of the relay nodes into 
two sets 14 and V£, such that \Vk\ = k, \V£\ = m — k, 
V k U Vfc c = K and V k n V k c = 0. Let 

h 

^=EE E E E c *> ^ 

then Cfc is the capacity of the corresponding S-t-cut, and 



E[C k ] = E x Ep[C k ] 

= [(m — k)h + fc(m 



fc) + fe]C, (47) 



Now, let Cs.t be the minimum cut capacity among all 
S-t-cuts, 



Cs,t = mm O fc . 

0<fe<m 



(48) 



By comparing © and d47l i, we note that we on longer 
have symmetry in E[Ck] with respect to k, i.e., E[Ck] ^ 
E[C m -k] for k = 0,1,..., m. In the single source case, 
the minimum value of E[Ck\, i.e., E[C s .t] is obtained when 
k = or k = to due to the symmetry (J5[Cq] = E[C m ]). 
This means that the bottlenecks are at the source end and also 
the destination end. Nevertheless, when we have multiple 
sources, £7 [Co] > E[C m ], and the minimum expectation 
value of the capacity among all cuts with any size is 
E\Cs,t] = E[C m ), which implies that we have only one 
bottleneck at the destination end. 

For the given set of source nodes S = {sx, s s } and the 
sets of destination nodes T = {ii, ...,£;} and relay nodes 
1Z = {ui, ...,u m }, define the network coding capacity for 
multiple sources and multiple destinations as 



Cs,r = minCs *. 



(49) 



Then, by the same method used in the previous section, we 
can show that Cs.t concentrates at E[C m ] = mC with high 
probability, where C is defined the same as before. This 
indicates that Cs t an d C s j concentrate at the same value. 
This is because they have one bottleneck in common. 

Theorem 11: When n is sufficiently large, the network 
coding capacity Cs,t satisfies 

Pr(Cs, T > (1 - e a )E[C m \) = 1-0 ( -L ) , (50) 

V to" / 
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Fig. 3. Interference at each node in G(X , Pn>7) 
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Fig. 5. Interference at each node in G(X , V, 7) 
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Fig. 4. Capacity of random s-t-cut of size k = 50 in G(X, Po,7) 
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Fig. 6. Capacity of random s-t-cut of size k = 50 in G(X, V , 7) 



where e a 
mC. 



_ (y+l)R 
- E[C m ] 



\/2am\xim for a > and E[C m ] = 



Proof: The proof is the same as that for Theorem [9] by 
replacing of E[C ] by E[C m ) □ 

Theorem 12: When n is sufficiently large, the network 
coding capacity Cs t satisfies 



Pr(C 5 . T < (l + e a )E[C m \) = 1-0 



1 



(51) 



where e a 
mC. 



^^jV2am\nm for a > and E\C m ) = 



Proof: The proof is the same as that for Theorem [10] by 
replacing of E[Cq] by E[C m ] □ 

V. Simulation Studies 

In this section, we present some simulation results on the 
SINR model and network coding capacity. Fig. 3 and Fig. 4 
show simulation results on interference and cut capacity in 
G(X,P ,i), where n = 2000, L{x) = ^x^, N = 
0.02, fi = 0.2 and 7 = 0.02, and every node transmits with 
constant power Pq = 0.01. Fig. 5 and Fig. 6 show simulation 
results on interference and cut capacity in G(X, V, 7), where 



n = 2000, L(x) = ^-x' 3 , N = 0.02, p = 0.2 and 
7 = 0.02, and every node transmits with power P uniformly 
randomly distributed over [0.01,0.02]. The results confirm 
the concentration behavior of interference and cut capacity. 

VI. Conclusions 

In this paper, we studied network coding capacity for 
random wireless networks with interference and noise. In 
this model, the capacities of links are not independent. By 
using coupling and martingale methods, we showed that 
when the size of the network is sufficiently large, the network 
coding capacity still exhibits a concentration behavior in 
cases of single source multiple destinations and multiple 
sources multiple destinations. We demonstrated simulation 
results that meet our theoretical bounds of network coding 
capacity. 
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